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FINITE-DIFFERENCE COMPUTATIONS OF ROTOR LOADS 
F. X. Caradonna and C. Tung 

Aeromechanics Laboratory, U.S. Army Research and Technology Laboratories — AVSCOM 

NASA Ames Research Center 
Moffett Field, California 


Abstract 

This paper demonstrates the current and 
future potential of finite-difference methods for 
solving real rotor problems which now rely largely 
on empiricism. The demonstration consists of a 
simple means of combining existing finite- 
difference, integral, and comprehensive loads 
codes to predict real transonic rotor flows. 

These computations are performed for hover and 
high-advance-ratio flight. Comparisons are made 
with experimental pressure data. 
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Notation 


x - x'/c, nondimensional streamwise coordinate 
(prime denotes a physical quantity) 

y « y'/R, nondimensional spanwise coordinate 
(prime denotes a physical quantity) 

z * z'/cfi 1 ^, scaled nondimensional normal 
coordinate (prime denotes a physical 
quantity) 

a - angle of attack, deg 

a p - partial angle of attack, deg 

Y * specific heat ratio 

5 - thickness ratio 


AR 

B 

C 

Cr 


R/c, aspect ratio 

2M 2 f/AR6 2/3 

1/AR 2 6 2/3 


lift coefficient 


9^ - collective pitch angle, deg 

0 1 - longitudinal cyclic input, deg 

0 1s - lateral cyclic input, deg 
U * advance ratio 


C T - thrust coefficient 
c ■ chord 
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potential 


a scaled velocity perturbation 


D - Bg 


$ * normal derivative of <J> 

n 
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1 - f M* 
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>x - m t 


f$ 

Y x x 


$ » second normal derivative of <j> 

nn 

i> * azimuthal angle, deg 


f - y + \i sin 


* angular speed 


g « x/AR - v cos ^ 


- solidity 


M t - tip Mach number 
R - blade radius 
r - radial station 
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Introduct ion 

The design of new rotor systems can be a 
risky undertaking. Therefore, most of these 
machines are derivatives or scaled versions of 
other well-proved rotor systems. Much of the risk 
of helicopter development stems from an inadequate 
understanding of and ability to predict the per- 
formance-def ining physical behavior of the rotor — 
especially the aerodynamic behavior. Most of the 
various flow phenomena that limit or define the 
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capabilities of a helicopter (e.g., wake forma- 
tion, dynamic stall, acoustic propagation, tran- 
sonic behavior, blade/wake interactions, and other 
unsteady loading and drag sources) are problems 
the understanding of which has not progressed far 
beyond the empirical level. The translation of 
our empirical understanding into a more reliable 
predictive ability probably awaits the development 
of more comprehensive numerical modeling tech- 
niques than currently exist. At present, our best 
computational schemes are built around various 
integral methods, of which lifting-line theory is 
the most commonly used. These methods are, in the 
main, limited to the prediction of linear compres- 
sible flows. The use of tabular airfoil data 
extends the validity of these methods, but only to 
the extent that two-dimensional, steady flow is 
representati ve of the rotor flow. 

A much more general predictive tool is to be 
found in the differential methods of flow model- 
ing. This class of schemes (which includes the 
finite-difference methods) requires a high invest- 
ment in computer capability but has no intrinsic 
physical modeling limitations. Finite-difference 
methods for the solution of various high-speed and 
high-lift problems are becoming commonplace in 
fixed-wing aircraft design. However, such methods 
have not yet been employed for rotary-wing work 
because of the greater flow complexity. The 
intent of this paper is to demonstrate a means 
whereby existing codes (both integral and differ- 
ential) may be combined to solve realistic rotor 
flow problems, which have heretofore been treated 
empirically. 


Choi c e of Flow Models for Var lo u s Flow Problems 


The helicopter's various flow problems can 
require a wide number of models. If one were to 
compute an entire vehicle with accurate modeling 
of all relevant phenomena nothing less than a 
massive Navier-Stokes (NS) solver would be 
required. But far more modest goals have con- 
founded us for years. The oldest of these is the 
rotor free-wake problem. All integral (yortex 
lattice, etc.) methods have failed to produce a 
reliably workable tool — and for unknown reasons. 

i o 

A few preliminary f Inite-dif f erence attacks ' on 
the problem have not done better, but this work is 
only in its infancy^ --The best current approaches, 
are stTlX’ those that use prescribed wake models 
Retreating-blade stall is another problem that is 
best handled empirically. The problem has been 
treated with NS codes, but these are not yet 
satisfactory because of inadequate turbulence 
modeling. One area which seems no longer to be 
beyond our grasp is the high-speed, advancing- 
blade problem. Separation does not usually occur 
here, and the essential physics is inviscid. 


Euler codes are appropriate for these flows, but 
may only be required when shocks occur that gen- 
erate appreciable vorticity. Flows with weaker 
shocks are much more efficiently treated by poten- 
tial methods. Even with this simplest of finite- 
difference models, it is not currently feasible to 
treat an entire rotor flow, and the problem 
becomes one of efficiently embedding a local 
finite-difference solution in a global solution 
obtained by other means (probably an integral 
method). Such a scheme will still be nonlinear, 
three dimensional, and unsteady and could greatly 
reduce or eliminate the need for airfoil tables. 


Potential Equations of Motion and Boundary 
Conditions 

The starting point for the various potential 
formulations is the mass conservation and 
Bernouli's equation, together with the assumption 
of irrotationali ty. This system of equations will 
conserve mass, energy, and entropy when solved in 
a divergence form^ that maintains the separate 
identities of density and potential. However, a 
number of significant nonconservative 1 * methods 
solve only for potential. Such schemes would only 
be conservative if combined with a shock-fitting 
scheme. (Note that conservation implies the abil- 
ity to accurately predict shock motion and the 
resulting load variations.) The simplest poten- 
tial equations involve various small-perturbation 
approximations. This class of equations all have 
the form 


A<P 




B« 




C* 


yy 


D* 


xy 


(1) 


where F is a nonlinear flux term that is a func- 
tion of This equation is derived by means of 

a small cross-flow assumption and only one nonlin- 
ear term (that contained in F) is retained. A 
fortunate consequence of the small-perturbation 
approximation is that the equation remains conser- 
vative in spite of it3 not explicitly involving 
density. 


The boundary conditions for Eq. (1) are shown 
in Fig. 1. The computation region shown is a box 
which contains the outer section of a rotor. The 
boundary conditions on the faces of this box are 
usually undisturbed flow conditions, except for 
the inboard face where two-dimensional flow is 
assumed. The surface tangency condition is 
expressed by specifying the surface slopes on a 
mean plane of the rotor. In order to represent 
vortices and vorticity sheets, interior planes of 
potential discontinuity are imposed. In Fig. 1, 
two such planes are shown. One plane, which is 
aft of the trailing edge, models the trailing and 
shed vorticity. Another sheet below the rotor 
models a vortex located at the edge of this sheet. 
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The model shown in Fig. 1 automatically takes 
account of the inflow produced by the vorticity in 
the computation region. However, the vorticity in 
the remainder of the total rotor flow field is not 
automatically accounted for. The simplest way to 
include this outer flow seems to be to compute a 
partial angle of attack, a p , which is a result of 
the vorticity in the outer region, and impose it 
on the mean surface boundary condition. In 
return, the loads computed by the finite- 
difference computation provide the basic vorticity 
distribution for the outer flow computation. This 
exchange of flow angle and load information 
between the inner and outer flow regions is the 
heart of a f inite-difference embedding scheme. 

This embedding will now be demonstrated for two 
different flight situations. 


Model Hover Testing 


The simplest problem that requires the above 
embedding procedure is a high-speed hover. 

This problem is relatively simple because it is 
steady, and the required validation data are 
fairly easy to obtain on the model scale. Fig- 
ure 2 shows a 7.5-ft-diam model which was tested 
in the Aeromechanics hover facility. The model 
was equipped with a hub-mounted Scanivalve assem- 
bly which permitted the measurement of pressure at 
120 points on the rotor. .Since the wake geometry 
is unpredictable/ its measurement is a vital ele- 
ment of such testing. These measurements were 
performed with traversing hot wires. The result- 
ing information was then used as the required 
prescribed wake data in an integral flow code — in 
this case the Analytic Method, Inc. code, HOVER. ^ 
For subcritical tip speeds, HOVER predicts the 
measured lift distribution very well. However, 
for supercr itical tip speeds, HOVER cannot treat 
the resulting transonic nonlinearity , and it 
becomes necessary to use an embedded finite- 
difference solution. For the present problem a 
steady rotor, small-perturbation code was used.® 
The boundary conditions employed were the same as 
depicted in Fig. 1, except that two vortices were 
included in the computation region. In the course 
of the computation, the locations of the vortices 
were fixed by the experimental data, and their 
strength was equated to the maximum blade bound 
circulation. A typical comparison of measured and 
computed pressure distribution is shown in Fig. 3 
for a tip speed M^. - 0.877. 


The fact that the f inite-difference embedded 
solution reproduces the measured loads well beyond 
the range of validity of the linear HOVER code 
demonstrates the validity of the small- 
perturbation method. It also shows that the 
scheme keeps consistent account of vorticity, in 


spite of this vorticity being split into two 
entirely separate codes. 


Forward Flight Computations 

Transonic flows are intrinsically unsteady, 
and advancing-blade transonic flows can be among 
the most unsteady flows of all. Nevertheless, 
tabular profile data take no account of this 
unsteadiness. Such flows are difficult to compute 
if the wake is to be accounted for, as was done in 
the previous hover case. For such cases, the 
vortices from previous blades are no longer neatly 
aligned with a coordinate; instead, they assume a 
wide range of orientations as they move through 
the grid. At present, such computations are being 
attempted only in two dimensions. However, a 
simple first approach to the problem is to account 
for the inflow entirely with a partial angle of 
attack, with no vortex in the grid other than the 
trailing and shed vorticity sheet. Such a scheme 
is almost certainly valid at high advance ratios 
where the induced inflow is small. This scheme 
has been implemented in a code 0 we call FDR 
(Finite-Difference Rotor Code); FDR solves the 
low-frequency variant of Eq. (1) (that is, A - 0). 

Another problem with forward flight computa- 
tions is that the blade is constantly responding 
to the varying loads, and a trim solution must be 
obtained. The problem here is one of efficiency. 
Trim solutions are iterative processes often 
requiring the rough equivalent of computing 
10 rotor revolutions. This much computing is no 
problem when the blade forces are obtained from a 
table, but it would be prohibitive for a finite- 
difference solution. A solution to this problem 
is to perform the basic f inite-dif f erence solution 
outside the central trim iteration. 10 Inside the 
trip loop, tables (or any fast, convenient approx- 
imate to a finite-difference solution) can be used 
to find a load correction based on the new blade 
motion. That is, 


C/a) 


l fd old 


♦ c. 


(a) - C 


table 


L ^old 5 
table 


where a and a Ql(i are the angles of attack from 
the current and previous trim loops, respectively. 
The solution is converged when a •+• and the 

lift correction vanishes; that is, when the 
finite-difference lift is fully consistent with 
the rotor inflow and motion. This process has 
been followed using the comprehensive rotor code 
CAMRAD 1 1 to provide the outer inflow, blade 
motion, and trim computations. 

The matching of the CAMRAD and FDR codes is 
summarized in the flow diagram shown in Fig. 4. 

The process is started by obtaining a trimmed, 
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nonuniforra inflow solution with the lift totally 
obtained from airfoil tables. This is, in fact, 
the normal operation of CAMRAD, except that par- 
tial wake-influence coefficients are computed in 
addition to the usual full-wake values. These 
partial coefficients are used to obtain partial 
angles of attack, which in turn are fed to the 
FDR. The lift distributions so obtained are then 
fed to the CAMRAD trim loop where the airfoil 
tables are used to find a lift correction. The 
process now iterates between the FDR and the trim 
loop until convergence of a is achieved. This 
scheme i3 an efficient one because it keeps all 
the time-consuming computations (the influence- 
coefficient and especially the finite-difference 
calculations) out of the innermost trim loop. It 
has been our experience to date that the conver- 
gence of this scheme is extremely rapid. 

A suitable body of experimental data on which 
to test the present code was obtained at a recent 
model-rotor test performed at the Deutsch- 
Niederlaendischer Windkanal (DNW). 12 The rotor is 
a 1 /7-scale model of the Cobra operational load 
survey (OLS) blades (Fig. 5). Like the OLS 
blades, the model is pressure instrumented (but 
not nearly so densely). At the present time, a 
small amount of the pressure data is available. 
Moreover, the model is fairly uncomplicated and is 
a good first test case for the present method. 

The rotor is a teetering rectangular blade with a 
linear twist and has a modified BHC-540 airfoil. 
Since the model is quite stiff, we have not con- 
sidered blade bending or twisting in the present 
computations. 


there is a moderate difference between the C# 

L FD 

and the C L . However, for N ■ 2 there is no 
difference between them, and the solution appears 
completely converged. Another view of the conver- 
gence is shown in Fig. 7, which shows the azi- 
muthal variation of a for various Iterations. 
Again the solution appears completely converged by 
the second iteration, with the greatest change 
occurring in the second quadrant. 

Figure 8 shows the spanwise variations of 
C L for a wide range of azimuths for the pure 
CAMRAD and the converged CAMRAD-FDR solutions. It 
is seen here that the effect of switching on the 
f inite-difference solution is to shift the loading 
somewhat from the rear to the front of the rotor 
disk. Accordingly, the greatest trim effect is on 
the lateral cyclic input (Fig. 9). Furthermore, 
the effect of the f inite-dif Terence solution is to 
bring the trim settings into somewhat closer 
accord with the experimental values. Note also 
that the trim settings are essentially fixed after 
the first iteration. 

At present, the only available pressure data 
are for the upper blade surface at r/R - 0.95. 

The comparison of the measured and computed upper- 
surface pressures is shown in Fig. 10. A similar 
pressure comparison at an advance ratio of 0.345 
is shown in Fig. 11. There is a considerable 
difference between measured and computed pressures 
in the first quadrant where the pressures appear 
to be generally underpredicted . The comparison 
tends to be better in the second quadrant, where 
the shock location is correctly predicted. 


The available data sets are high-speed cases. 
A typical run, which we have computed, has a tip 
Mach number My of 0.663# an advance ratio y of 
0.298, and a C^/o - 0.0769. For this case, the 
contribution of the grid-enclosed trailing vortex 
to the full angle of attack is seen in the span- 
wise distributions of a and a at i|> - 90° 

(Fig. 6). At this azimuth, a Eip vortex passes 
beneath the rotor at 0.7R. Accordingly, a is 
augmented outboard and diminished inboard of this 
location. Since the rotor lift varies similarly, 
the trailing vorticity will have an opposing 
effect on the induced flow of the vortex. Remov- 
ing this opposing effect results in an a dis- 
tribution with a much stronger variation. P This ( 
large difference between a and a demonstrates 
the importance of keeping a consistent account of 
the wake with a matched scheme. In Fig. 6, the 
distribution of a, a 


L FD 


(lift computed in 


FDR), and (the lift computed in the CAMRAD 
trim loop) are shown for three iterations. The 
N - 0 iteration is the starting CAMRAD lifting- 
line solution. All subsequent iterations are 
complete integral-differential cycles. For N - i 


Concluding Remarks 


The above results are an indication of the 
most probable immediate direction of finite- 
difference computations for application to high- 
speed rotor loads. The mixed integral- 
differential approach to transonic load prediction 
has the advantage of making immediate use of our 
existing loads codes. It also can greatly reduce 
the need for two-dimensional airfoil testing. (It 
should be borne in mind that two-dimensional pro- 
f lie data are proba b 1 y o f dubfous^val ue~ when non-" 
rectangular planforms or t r a ns o n i c" uns t e adl n es s 
are'irnvolved.y^^n^ts^resent stage of develop- 
ment, the approach gives fair to good comparison 
with experimental pressure data. However, consid- 
erably more comparisons with forward flight data 
are required. The method clearly can be enhanced 
by employing a conservative, full-potential 
f inite-difference model and by improving the near- 
wake modeling capability. This latter improvement 
would involve inserting the moving vortex back 
into the finite-difference grid. Such 
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modifications will certainly be required if the 
method is to be extended to lower advance ratios 
at which the blade/vortex interactions are more 
important .Even in its present form, however, the 
mixed integral-dif ferential scheme affords a 
considerable improvement in our flow prediction 
capabilities. 

In the future, f inite-dif ference modeling can 
be developed and extended for application to a 
wider variety of problems, the most important of 
which is probably the prediction of rotor-wake 

behavipr.. Ultimately, it-wilJL not be unreasonable 

to expect the development of unified finite- 
difference models of the entire rotor/fuselage. 
However, this goal will require many fundamental 
advances in numerical modeling and in basic flow 
physics. „ 
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FAR FIELD BOUNDARIES 



Fig. 1 Boundary conditions for the rotor potential-flow problem. 
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MEASUREMENT 



a) r/R - 0.89 and 0.96. 

Fig. 3 Comparison of measured and computed chordwise pressure distributions for hover model: 

- 0.877, AR - 6.0, 9 c * 8°. 
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M t - 0.877, AR - 6.0, 6 e - 8° 
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b) r/R - 0.68 and 0.80. 
Fig. 3 Concluded. 
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Fig. 4 A hybrid integral-differential rotor flow scheme. 
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Experimental setup of the model rotor at DNW 










(a) N » 0 


<b) N » 1 


(c) N * 2 


Fig. 6 Evolution of the load and angle-of -attack distribution with CAMRAD-FDR iteration number 

M t - 0.663, u - 0.298, ty - 90°. 


TIP MACH NO. - 0.663 
ADVANCE RATIO - 0.298 



Fig. 7 Azimuthal inflow variation as a function of iteration number: Mm » 0.663, u * 0.298. 
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Fig. 10 Comparison of computed and experimental surface pressure for an advancing rotor: ftp 

u - 0.298, r/R - 0.95. 
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